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We investigate eigenvalues of many-body systems interacting by two-body forces as well as those 
of random matrices. We find a strong linear correlation between eigenvalues and diagonal matrix 
elements if both of them are sorted from the smaller values to larger ones. By using this linear 
correlation we are able to predict reasonably all eigenvalues of given shell model Hamiltonian without 
complicated iterations. 
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Diagonalization of matrices is a very common prac- 
tice in many fields. There have been a number of algo- 
rithms for diagonalization, such as the Jacobi method, 
the Householder method, the Lanczos method, etc. 
Full diagonahzations of obtaining all the eigenvalues be- 
come difficult virhen the dimension of matrices is larger 
than 10^. One can obtain a very few lowest eigenvalues 
of sparse matrices when dimension of matrices is 10^"^*^ 
by super computers. There have been also many efforts 
in obtaining eigenvalues by using the energy centroid and 
spectral moments, and typical works along this line can 
be found in, e.g., Refs. H, i, i, i, i]. Unfortunately, 
the eigenvalues such obtained by using centroids and mo- 
ments are valid at the statistical level (i.e., after the en- 
semble average) and not good enough for individual di- 
agonahzations. 

The purpose of this Letter is to report a remarkable 
correlation between eigenvalues and diagonal matrix el- 
ements, if both are sorted from smaller values to larger 
ones. One can predict approximately all eigenvalues of 
given matrices by using this correlation. Although eigen- 
values such predicted are not exact, they are good ap- 
proximations of exact results (within a precision of 1 — 4% 
in examples that we have studied). 

Let us begin with Hamiltonian matrices which con- 
serve the angular momentum. The two-body interac- 
tion parameters are taken to be Gaussian random num- 
bers (i.e., two-body random ensemble which is called 
the TBRE). In such cases both the diagonal matrix ele- 
ments of a many-body system and their eigenvalues are 
known to exhibit Gaussian distributions approximately. 
One easily conjectures that there might exist a linear 
correlation between these two sets of quantities. Let 
us denote exact eigenvalues obtained by diagonalization 
by using E^'^'^'^^, dimension of matrix by D (number of 



spin / states is denoted by Dj), predicted eigenvalues 
by i<;exact^ diagonal matrix elements by Hu, and aver- 
age energy by H. We define the linear correlation co- 

ethcient r = — _ „ _ ■ the absolute 

value of r, |r|, is less than or equals to 1. If there ex- 
ists a strong linear correlation between Ef^^'^^ and Hu 
{i — 1,2, ■ ■ -n), \r\ — > 1. We investigate below whether 
or not the coefficient r is very close to 1 by a number 
of different configurations. We take six successive sets 
of random interactions for each system in demonstrating 
the correlation between eigenvalues and diagonal matrix 
elements. 

In Fig. l(a-f), we present exact eigenvalues versus 
diagonal matrix elements of J = 20 states (dimension 
Dj ~ 29) for four identical fermions in a single-j {j — ^) 
shell. Both eigenvalues and diagonal matrix elements are 
sorted from the smallest to the largest. The linear corre- 
lation coefficients r are more than 0.95 in all these panels. 
Each panel corresponds to one set of random interaction 
parameters. 

In Fig. 2(a-f), similar results are presented for sd bo- 
son systems [8(]. Here 1 = 6 and boson number n = 12 
(dimension Di = 37). In Fig. 3(a-f) we present results 
for the I = 17 {Dj ~ 508) states of three valence protons 
and three valence neutrons in the (2si/2j 1^3/27 0«ii/2) 
shell; in Fig. 4 we present results of / = 4 {Dj = 3017) 
states of the same configuration. The Hamiltonian that 
we use in Figs. 3-4 conserves parity and isospin as well 
as angular momentum, i.e., here we use the conventional 
shell model Hamiltonian in nuclear structure physics. 
The two-body interactions are taken to be random values 
following the Gaussian distribution. Correlations in all 
these figures are seen to be striking, although there are 
small deviations when eigenvalues are around the small- 
est or the largest. 

Now we discuss how one obtains the predicted linear 
correlation by a simple algorithm. Let us assume that 
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where En is the eigenvalues of the Hamiltonian. By using 

D D 



we obtain 
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FIG. 2. (Color online) The same as Fig. 1 except that 
these examples correspond to sd-boson systems. Here 
boson number n = 12, J = 6, Dj = 37. 
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B = {1- A)H. 
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All straight lines in this paper are plotted by using the 
above relation. 
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FIG. 1. (Color online) Linear correlation between 
exact eigenvalues and diagonal matrix elements. These 
examples correspond to four fermions in a single-j shell 
with j = 31/2. Here / = 20, Di^29. r is the linear 
correlation coefficient. Straight lines are plotted by using 
Eq. (1), with A and B values evaluated by Eq. (3). 
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FIG. 3. (Color online) The same as Fig. 1 except that 
these examples correspond to three valence protons and 
three valence neutrons in the (2si/2, 1^3/2, 0«ii/2) shell. 
Here / = 17 and £>/ = 508. 

In the above examples the two-body interacting pa- 
rameters are given by Gaussian random numbers. The 
matrices of these systems are given by coefficients of frac- 
tional parentage determined by geometry of the configu- 
ration and two-body interactions which are taken to be 
Gaussian. Now we come to matrices in which all matrix 
elements are given by random numbers. 

In the first set of examples, all matrix elements are 
taken to be uniformly distributed random numbers be- 
tween — 1 and 1. Figure 4 shows exact eigenvalues ver- 
sus diagonal matrix elements (dimension £'=500), both 
sorted from smaller to larger values. The correlation is 
seen to be very good, except for small deviations near 
the largest and the smallest eigenvalues. 
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FIG. 4 (color online) Eigenvalues versus diago- 
nal matrix elements {di — 500). Both eigenvalues and 
diagonal matrix elements are sorted from smaller values 
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to larger ones. All matrix elements are taken to be uni- 
formly distributed random numbers. One sees that the 
eigenvalues predicted by linear correlation and the exact 
values are almost equal except for a very few around the 
maximum and minimum values. 

We also investigate the emergence of the above corre- 
lation in matrices for which all matrix elements are ran- 
dom numbers following the Gaussian distribution. Fig- 
ure 5 shows results of matrices with D = 500 and all 
matrix elements being Gaussian random numbers. The 
correlation seems not linear here, but rather a hyperbolic 
tangent function because eigenvalues in these case follow 
Wigner's semicircle distribution. However, it is impor- 
tant to point out two important facts. (I) There also ex- 
ists a remarkable linear correlation between eigenvalues 
and diagonal matrix elements on average in these cases. 
The above linear correlation is actually very good except 
for eigenvalues near the smallest or largest ones, because 
the level density of eigenvalues in the middle part is much 
larger than that near the minimum or maximum. All 
"medium" eigenvalues can be predicted very well by the 
linear correlation determined by Eqs. (1-3). For matrices 
with dimension larger than 20, correlation coefhcient r is 
larger than 0.95. We note without details that one eas- 
ily obtains the hyperbolic tangent function which seems 
proper in fitting the results around the minimum or max- 
imum eigenvalues, by using matrix elements in the same 
way as we do for the TBRE systems. (II) The correla- 
tion coefficients between eigenvalues and diagonal matrix 
elements are larger than 0.995 if one artificially enlarges 
diagonal matrix elements by about two orders or more 
(typically by a factor of 50-60) for Gaussian random ma- 
trices. The latter fact is very important, because in the 
nuclear shell model calculations, magnitudes of diagonal 
matrix elements are usually much larger than (non-zero) 
off-diagonal matrix elements by about one or two orders. 



much larger than that of eigenvalues near the minimum 
or maximum. All r values are larger than 0.98. 

Now we exemplify our approach by a realistic example. 
In Fig. 6(a), we present eigenvalues of spin / = states 
of ^"'Mg {Dj — 1161) versus the diagonal matrix elements 
of the same configurations. Here we use the USD inter- 
actions It is seen that the linear correlation is very 
good. Similar results can be found for other spin states. 
In Fig. 6(b) we show exact eigenvalues Ef^'^'^^ in compar- 
ison with those predicted by using linear correlations {Ei 
in Eq.(l)). Although eigenvalues predicted by using the 
linear correlation deviates from the exact values slightly, 
they are close to the exact ones in general. More impor- 
tant, the agreement between predicted eigenvalues and 
exact results is substantially improved when one goes to 
the medium region. As shown in Fig. 6 (c), predicted 
100th-105th 0+, 2+, 3+, 4+, 6+ states are very close to 
exact eigenvalues obtained by diagonalizations. 
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FIG. 5. (Color online) Correlation between eigen- 
values of random matrices (Gaussian type) and diagonal 
matrix elements. The dimension is 500 here. It is noted 
that the linear correlation is actually very good except 
for eigenvalues near the smallest or largest ones, because 
the level density of eigenvalues in the middle region is 



FIG. 6 Energy Spectrum for '^^Mg. The USD inter- 
actions are used here, (a) Eigenvalues of 7 = states 
obtained by diagonalization versus diagonal matrix ele- 
ments. One sees a striking linear correlation. The cor- 
relation coefficient r is larger than 0.999. (b) Energy 
spectrum by exact diagonalization and that predicted 
by linear correlation. The eigenvalues of the lowest two 
1 — 2 levels (and two lowest 7 = 6 levels) are degener- 
ate because diagonal matrix elements corresponding to 
these states are equal due to the isospin conservation of 
the USD interactions, (c) Energies of excited states in 
medium region. The left hand side corresponds to ex- 
act energies obtained by diagonalizations, and the right 
hand side corresponds to predicted values by using linear 
correlations. The 100th-105th 0+, 2+, 3+, 4+, 6+ states 
are presented. 

To summarize, in this paper we have discovered the 
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linear correlation between exact eigenvalues and diagonal 
matrix elements for individual runs of the two-body ran- 
dom ensemble or random matrices, if both exact eigen- 
values and diagonal matrix elements are sorted from the 
smaller values to larger ones. We exemplify this correla- 
tion by fermions in a single-j shell, many-j shells and sd 
bosons by using random interactions, and energy spectra 
of ^''Mg by using the USD interactions. 

This linear correlation seems to be universal. For two- 
body random ensembles, the linear correlation coefficient 
is ~0. 95- 1.00; for matrices in which all matrix elements 
are uniformly distributed random numbers, the linear 
correlation is very good; for Gaussian type random ma- 
trices, there are sizable deviations between exact eigen- 
values and predicted eigenvalues only if eigenvalues are 
close to the minimum or maximum, and in this case the 
linear correlation is well applicable to predicting levels in 



the medium region. 

Although eigenvalues obtained here are not exact, they 
are very good approximations of exact values. We there- 
fore believe our results are very useful in studying prob- 
lems in which one needs approximate eigenvalues when 
other methods are not possible. 
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